Correlated electron transport through double quantum dots coupled to normal and 

superconducting leads 
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We study Andreev transport through double quantum dots connected in series normal and super- 
conducting (SC) leads, using the numerical renormalization group. The ground state of this system 
shows a crossover between a local Cooper-pairing singlet state and a Kondo singlet state, which is 
caused by the competition between the Coulomb interaction and the SC proximity. We show that 
the ground-state properties reflect this crossover especially for small values of the inter-dot coupling 
t, while in the opposite case, for large t, another singlet with an inter-dot character becomes domi- 
nant. We find that the conductance for the local SC singlet state has a peak with the unitary- limit 
value 4e 2 /h. In contrast, the Andreev reflection is suppressed in the Kondo regime by the Coulomb 
interaction. Furthermore, the conductance has two successive peaks in the transient region of the 
crossover. It is further elucidated that the gate voltage gives a different variation into the crossover. 
Specifically, as the energy level of the dot that is coupled to the normal lead varies, the Kondo 
screening cloud is deformed to a long-range singlet bond. 

PACS numbers: 73.63.Kv, 74.45.+C, 72.15.Qm 



I. INTRODUCTION 

A quantum dot (QD)i is a prototype of nanoscale sys- 
tems, and has offered unprecedented opportunities to un- 
cover correlation effects on quantum transport. In par- 
ticular, the Kondo effect in a QD system gives rise to re- 
markable properties in electron transport^— which has 
encouraged further studies in this field. With recent ex- 
perimental advancements, it becomes possible to exam- 
ine the Kondo physics in a variety of systems, such as an 
Aharonov-Bohm (AB) ring with a single QD and double 
quantum dots (DQD). 

A QD coupled to superconducting leads is also an 
intriguing system, which can bring about the competi- 
tion between superconductivity and the Kondo effect. 
In fact, such a competition has recently been observed 
in carbon nanotube and semiconductor QD systems.— ~— 
Furthermore the Andreev scattering at a junction of 
a normal metal (N) and superconducting (SC) leads, 
is also a fascinating phenomenon. Particularly, in the 
QD's coupled to the normal and SC leads a short-range 
Cooper pair penetrates into the QD as a linear com- 
bination of the empty and doubly occupied states, so 
that the SC proximity becomes sensitive to the Coulomb 
repulsion. The interplay between the Andreev scatter- 
ing and the Kondo effect takes place in such a situa- 
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FIG. 1: (Color online) Quantum dots (QD's) coupled to 
normal (N) and superconducting (SC) leads in series. 



tion, and has been studied intensively for a single QD, 
theoretically^— and experimentally . 25 ' 26 We have stud- 
ied the ground-state properties of a single dot Andreev- 
Kondo system, and have confirmed in a previous work 2 ^ 
that the zero-temperature conductance shows a maxi- 
mum at a crossover region between the Kondo singlet 
state and the local Cooper-pairing singlet state. The lo- 
cal Cooper pairing is caused by the SC proximity effect, 
and consists of a linear combination of the empty and 
doubly occupied states. We shall refer to this singlet 
state as the local SC singlet state for short. 

The Andreev transport has also been explored further 
in other nanosystems, e.g., an AB ring with a single 
QD;22 a DQD system^ and molecular wires32 Further- 
more, the Josephson current in a serial DQD system, 
where both leads connected to dots are in SC states, 
has also been investigated in recent years i 30 ' 31 It is clar- 
ified that the inter-dot coupling plays a crucial role on 
the Josephson effect. In fact, one of the typical features 
of the DQD systems is that another type of a singlet 
ground state, which is referred to as an inter-dot singlet 
hereafter, appears and plays an important role on the 
low-energy properties. To our knowledge, however, there 
have been few studies of the Andreev transport through 
the DQD system so far, while a lot of studies have been 
done for the normal transport in the DQD system. The 
underlying Andreev-Kondo physics in the DQD system 
is still less understood, and is needed to be clarified pre- 
cisely. Moreover, the DQD system coupled to normal and 
SC leads can be fabricated ) 32 ' 33 and is studied from the 
viewpoint of spin entanglement^ These developments 
are also spurring further research of the DQD system 
coupled to normal and SC leads. 

In the present paper, we theoretically study the trans- 
port through the DQD coupled in series to normal and 
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SC leads, as shown in Fig. [TJ The SC proximity into 
the DQD is affected more by the Coulomb interaction L/ 2 
at the QD2, which is located adjacent to the SC lead, 
than the interaction U\ at the QD1, sitting away from 
the SC lead (see Fig. [1]). Such a spatial variation gives a 
variety, which cannot be realized in the single dot, into 
the Andreev-Kondo physics in the DQD system. Par- 
ticularly, we focus on the crossover between the ground 
states which can be classified into the local SC singlet, the 
Kondo singlet and the inter-dot singlet states, and clar- 
ify how the crossover affects the transport properties. We 
first of all demonstrate that in the electron-hole symmet- 
ric case this system can be mapped onto the two impu- 
rity Anderson model coupled to normal leads. Then, we 
calculate the conductance due to the Andreev reflection 
in a wide parameter range of the inter-dot coupling, the 
Coulomb interaction, and the energy level of the quan- 
tized states in the dots, using the numerical renormaliza- 
tion group (NRG) method^^ in the limit of a large SC 
gap A — > oo. In the case that the Coulomb interaction 
XJ<x at the QD2 is small, the short-range Cooper pair can 
penetrate into the QD2, and the conductance takes the 
unitary-limit value 4e 2 //i. In contrast, for large [/ 2 the 
SC proximity to the QD2 is suppressed, and the conduc- 
tance docs not reach the unitary-limit value. It reflects a 
significant change of the SC correlation penetrating into 
the DQD, which we have deduced from the behavior of 
the renormalized parameters for the Bogoliubov parti- 
cles. Furthermore, we find that the conductance has two 
peaks with the unitary-limit value near the crossover be- 
tween the Kondo singlet state and the local SC singlet 
state. The contrast between the two singlet states is pro- 
nounced for small inter-dot coupling t. In the opposite 
limit, for large t, the ground state becomes another sin- 
glet with an inter-dot character. We examine also the 
gate voltage dependence of the conductance, varying the 
energy level of each dot separately. The crossover be- 
tween these singlet states takes place as the energy levels 
are varied, and near the crossover point the conductance 
shows a peak. It is further found that a different type 
of Kondo singlet with a long-range singlet bond emerges, 
as the energy level of the QDI which is located at the 
normal-lead side varies. 



This paper is organized as follows. In the next section, 
we introduce the model and give a brief outline of the 
Bogoliubov transformation. Then in Sec. IIII1 we show 
the numerical results and discuss the Andreev transport, 
focusing on the interplay among the Kondo effect, the SC 
correlation and the inter-dot coupling. Moreover, the re- 
sults for the energy-level dependence of the conductance 
are presented. A brief summary is given in the last sec- 
tion. 



II. MODEL AND FORMULATION 
A. Model 

The Hamiltonian of a DQD coupled in series to normal 
(N) and superconducting (SC) leads is given by 

H = Hdqd + H$ + Hn + Ht,s + Ht,n, (1) 

where Hdqd arid Hs(n) represent the DQD part and 
the SC (normal) lead part, respectively. Ht,s(n) is the 
mixing term between the QD and the SC (normal) lead. 
The explicit form of each part reads 

Hdqd = ^ {& (n d)i - 1) + — ' (n d ,i - l f\ 
t=l,2 ^ ' 

a 

H S =J2 £ S,kcl M c SM - ( A 4,fct 4,-/4 + H - C ') ' 
k,a k 

= ^ £ N,kC^ N _ ka C N ^ kal 
k.a 

H T ,N = Y.~jk { C N,kJlv + A* c Nm) > 

fc,<T 

Ht >8 = E % (osm*** + d lo SM ) • (2) 

k.a 

Here, £j = £j + t/j/2 for i = 1,2. The operator d[^ a 
creates an electron with energy £1(2) and spin a at the 
QD1(QD2). C/i(2) is the Coulomb interaction, and t is 

the inter-dot coupling, and n d ,i = J2 a d L d ia- c s(N),ka 
denotes the creation operator of an electron with the en- 
ergy £s(jv) ; fc in the SC (normal) lead. Vs/n is the tun- 
neling matrix element between the QD2/QD1 and the 
SC/normal lead, and A is an s-wave BCS gap. We as- 
sume that T s / N (e) = nVg, N J2k 5{e—Sk)/N is a constant 
independent of the energy e, where Af is the number of 
the states in each lead. 

To be specific, we concentrate on a large SC gap limit 
A — > 00 in the present paper. In this limit the quasi- 
particlc excitations in the continuum energy region above 
the SC gap are projected out. Nevertheless, the An- 
dreev reflection takes place inside the SC gap, and the 
essential physics of the low-energy transport is preserved 
still in the large gap limit.— In this case, the starting 
Hamiltonian H can be mapped exactly onto a single- 
channel model for which the NRG approach works more 
efficiently^^ 

H cS =Hg S + Hdqd + Hn + Ht.n , (3) 

Hf = -A d2 (4 f 4 + H.c.) , (4) 

A d2 = T s . (5) 

Here, an additional term Hg S appears instead of the SC 
lead Hs + Ht,s- The SC proximity effect becomes static 
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in the large gap limit, and is described by the pair poten- 
tial in the QD2, = ^s- This term breaks the charge 
conservation of the electrons, and causes the Andreev re- 
flection. 



B. Bogoliubov transformation in the £1 = case 

The system described by H has a conserved charge 
in the case that the value of the energy level of the QD1 
satisfies the condition ^ = 0. Then the Hamiltonian H 
can be transformed into an asymmetric two-impurity An- 
derson model for the Bogoliubov particles, the total num- 
ber of which is conserved , 24 i 28 1 40 This is due to the ro- 
tational symmetry in the Nambu pseudo-spin space, in 
which and A^2 can be regarded, respectively, as the z 
and x components of an external field fj which couples 
to the pseudo spin.— Specifically for £i = 0, the external 
field becomes finite only in the QD2, and its contribution 
to the energy is given by fj ■ f = — A^t^ + £,2 T z, where 
Tj for j = x,y,z is the Pauli matrix in the Nambu rep- 
resentation. Thus the system has a uniaxial symmetry 
along the direction of the local external field fj in the 
pseudo-spin space. 

In order to use these symmetry properties, we rewrite 
the normal lead part of the Hamiltonian, Hm and Ht.n, 
in a tight-binding form 



Here, / 0(T = J^ fc c N ka /V7\f. Note that no approximation 
has been made to obtain Eq. (j6]), and the hopping matrix 
element can be generated from EN,k via the House- 
holder transformation^ The tight-binding form reveals 
the electron-hole symmetry of Hm explicitly. In the case 
£1 = 0, the Hamiltonian can be simplified by the Bogoli- 
ubov transform, taking the direction of the local exter- 
nal field at the QD2 as a new quantization axis in the 
pseudo-spin space, 
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Hn = J2Y2 t N,n [fl+lafna + fLfn+la) > ( 6 ) 
n— cr 

H t ,n = £ V N ( fl d Xa + d\ a f Qa ) . (7) 
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FIG. 2: (Color online) Schematic of the energy level profile 
at the DQD part for the system after the Bogoliubov trans- 
formation. The (red) dashed line denotes the Fermi energy 
of the leads (ef = 0), and the (blue) dotted line indicates 
E2, which is located at the center between two energy levels, 
E2-U2/2 and£ 2 + £/ 2 /2. 



and ijv,-2 = t. Then the effective Hamiltonian H can 
be expressed in the form 



tf cff = E 2 (n 7 ,_ 2 - 1) + E f - 

2=1,2 

00 

+ EE l N,n (ll+lalna + H.C.) . (11) 



n=-2 a 

Here, n 7 ._; = J_ ia j_ i(7 is the number of Bogoliubov 
particles at the QDi for i = 1,2. In this representa- 
tion, the value of E 2 — U2/2 corresponds to an energy 
level for the Bogoliubov particles in the QD2. This cor- 
respondence is also illustrated schematically in Fig. [5] 
For instance in the atomic limit t N n — > 0, the QD2 
tends to be occupied by a single Bogoliubov particle for 
E2 — U2/2 < 0, while the QD2 tends to be empty for 
E2 — U2/2 > 0. The equation (TiTj) clearly shows that the 
total number of Bogoliubov particles Qtot is conserved, 



Qtot = 



00 

E 5 

J=-2 



'1,3 



(12) 



It should also be noted that the Friedel sum rule holds for 
the number of the Bogoliubov particles in the DQD ) 41 ' 42 



i=l,2 



2ip 

7T 



(13) 



Here, <p is the phase shift of the Bogoliubov particles, 
defined in the appendix [A] The phase shift ip can be 
deduced from the low-lying eigenvalues of discretized 
Hamiltonian of the NRG at the fixed point 4^ 

Furthermore, one can deduce the following quantities 
from the value of Q, via the inverse transformation of Eq. 



E> 



A 2 

^d2 



(10) 



The transformation is carried out for the whole cites, n > 
—2, including the DQD part for which we use a notation 
f—ia = di a for i = 1, 2. Similarly, the hopping matrix 
elements ijv n for n < are defined to be tpf _i = Vn 



A 



K= ^(-lfM =-- f(Q-2), (14) 

i=X,2 2 

M = ECW- 1 ) =^(Q-2). (15) 
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Here, Ka,i = d^d^ 



d li d lV 



The SC correlation K and 



the number of the original electrons M in the DQD cor- 
respond, respectively, to the average value of the x and 
z components of the induced pseudo-spin moment. 

In the case of £1 = 0, the linear conductance G at 
zero temperature can also be expressed in terms of the 
local charge Q, or the phase shift if, of the Bogoliubov 
particles (see also the appendix [Al . 



C. Conductance for /0 

For / 0, however, the uniaxial symmetry in the 
pseudo-spin space is broken, and the system no longer 
has a conserved charge such as the total number of the 
Bogoliubov particles Q to t ■ We calculate the conductance 
for these cases, using the Kubo formula^ 



G = lim ^ — \{GS\J N \n)\ 2 8(huj ~ e n ) 



Jn 



ie ^ V N / t 



N.kcr 



-JV,fco-"l<T 
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FIG. 3: (Color online) (a) Conductance G and (b) the pair 
correlations K as a function of t/Fjv for Ui = Uz = U and 
£i = £2 = 0. The hybridization energy scale is chosen to be 
Fjv = Ts, and the SC proximity appears through A<j2 = Ts- 
The NRG calculations have been carried out for A = 3.0 and 
T N /D = 1.0 x 10" 3 . 



Here, e n and |n) denote an eigenvalue and the corre- 
sponding eigenstate of H cS , |GS) is the ground state, 
and Jn is the current which flows from the normal lead 
to the QD1. 



III. NUMERICAL RESULTS 

In this section, we provide the NRG results for the 
ground-state properties of the serial DQD coupled to nor- 
mal and SC leads. The calculations have been carried 
out, taking the tunneling matrix elements tn,n for n > 
and V}v in the forms 



D 



1 + 1/A 



1 - 1/A 



n+l 



v N = 



y/1 - 1/A 2 »+Vl - 1/A 2 "+ 3 ' 

1 /1 + 1/A 



(19) 



2T N DA A 



2 \ 1 - 1/A 



log A, (20) 



where D is the half- width of the conduction band. We 
have kept lowest 1000 states in each NRG step, and have 
chosen the discretization parameter to be A = 3.0. 



A. Dependence on the inter-dot coupling t 

We first consider the case U\ = U2 = U, and choose the 
values of energy levels to be £ 1 =£2 = 0. In this case the 
energy level for the Bogoliubov particles in the QD2 is 
given by E- 2 = (= Tg), and the conductance is deter- 
mined by the local charge Q through Eq. (fTo| . The NRG 
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FIG. 4: (Color online) Conductance G as a function of i/IV 
near U /Ts = 2.0, where the other parameters are the same as 
in Fig. [3] This figure corresponds to the enclosed area drawn 
by the dotted line in Fig. 



results of the conductance are shown in Fig. [3] as a func- 
tion of the inter-dot coupling t for several values of U, 
for a fixed hybridization strength Tjv = ^s- For weak in- 
teractions U/Ts < 2.0, the conductance peak reaches the 
unitary-limit value 4e 2 //i, and the peak position shifts to- 
wards the small t side. In contrast, for U /Tg > 2.0, the 
conductance peak does not reach the unitary-limit value. 
Note that the onsite level for the Bogoliubov particle in 
the present case is given by E2 — U2/2 = Tg — J7 2 /2. Thus, 
for Ts < U/2, the Kondo effect takes place at the singly 
occupied QD2, as we will check out further in the next 
section. We can see that there is a fine structure near 
U /Ts = 2.0 for small values of t. Figure 2] is an enlarged 
picture of the conductance for the small t region, which 
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is marked with the dotted line in Fig. [3fa) . The con- 
ductance for U/Ts = 2.01 has two unitary-limit peaks. 
These two peaks get close to each other for U /Ts = 2.02, 
and the conductance shows a plateau for U/Ts = 2.04 
at 0.16 < t/Tfq < 0.3. Then, as U increases further, the 
maximum value of the conductance does not reach the 
unitary-limit value. 

In the present case, £i = £2 = 0, the SC pair correla- 
tion K defined in Eq. fTJJ and the conductance G in Eq. 
(|16p can be expressed in the simplified forms, 



K = -(it - (p), 

IT 



4 P 2 

G = —sin 2 2tp . (21) 
h 



Furthermore, the number of the real electrons takes a 
value of M — at half-filling. Note that the conductance 
takes the unitary-limit value 4e 2 //i for K = 1/2, which 
corresponds to ip = 37r/4 and Q = 3/2. In Fig. \3[h), the 
NRG results of K are plotted as a function of t/T^ for 
several values of U/Ts- We can see that the behavior of 
the pair correlation K also changes significantly at the 
Coulomb interaction of the value of U — 2T$- In the 
large t limit, two Bogoliubov particles occupy a bonding 
orbital consisting of the QD1 and QD2, so that Q —> 
2.0 and K — > 0.0 at t/T^ 1.0. For weak interaction 
U < 2Ts, K increases monotonically as t decreases, and 
approaches the value of K = 1.0, which corresponds to 
Q = 1.0, in the limit of t — >• 0. In contrast, for large 
interactions U > 2Ts, K shows a peak at an intermediate 
value of t, and then tends to K = 0.0 as t decreases. The 
conductance has two peaks in the case that the maximum 
value of K becomes K > 1/2. The behavior at U < 2Ts 
is determined by the Andreev scattering due to the finite 
SC pair potential Ad2 = Ts, while for large U > 2T$ the 
behavior is determined by the Kondo effect. 

It seems to be meaningful to compare these results with 
the normal transport through the serial DQD coupled to 
two normal leads (N-DQD-N) j£r£& particularly with the 
NRG results by Izumida et ai4& In the normal case, the 
conductance peak shifts towards the small t side as the 
Coulomb interaction U increases, and this behavior is 
similar to the one we observed for the Andreev trans- 
port in Fig. [3^a). However, in the N-DQD-N case, the 
peak height of the conductance reaches the unitary-limit 
value 2e 2 /h even in the large U region. Therefore, the 
suppression of the conductance for large U and the two- 
peak structure observed for small t, as mentioned in the 
above, are typical of the Andreev transport. 



B. Effects of the Coulomb interaction in each dot 

In order to study the difference between the Andreev 
behavior at U < 2Ts and the Kondo behavior at U > 
2Ts, we next consider the case that the Coulomb inter- 
action is switched on only in one of the two dots. In this 
situation, the role of the correlation in each dot can be 
seen separately. 
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FIG. 5: (Color online) Conductance for (a) Ui = and (b) 
U2 — as a function of the inter-dot coupling t. We set 
£i=6 = and F N = T s . 



In Fig.[5]the conductance for (a) U± = and (b) U2 = 
is plotted as a function of t/T^ for several values of (a) U2 
and (b) Ui, respectively. We can see in Fig.[5ja) that the 
conductance peak shifts towards the small t region as U2 
increases from to 2Fs. Then, for U2 > 2. Org the con- 
ductance decreases, and the peak height becomes smaller 
than the unitary-limit value 4e 2 /h. These features are 
similar to those we saw in Fig. for Ui — U2- In con- 
trast we can see in Fig.[5](b) that the peak position shifts 
towards the large t region as U\ increases keeping the 
peak height Ae 2 /h unchanged. Since U2 — in this case, 
the SC correlation can penetrate into the QD2 without 
being disturbed by the Coulomb interaction. Therefore, 
the QD2 can be regarded effectively as a part of the SC 
host. The inter-dot coupling t transmits the SC proxim- 
ity from the QD2 to the QD1. For this reason, even in 
the presence of the Coulomb interaction Ui, the conduc- 
tance can take a maximum with the unitary-limit value 
4e 2 //i. It explains the reason why the conductance peak 
shifts towards the larger value of t as U\ increases. We 
see also in Fig. EJb) that there are no pronounced qual- 
itative changes in the feature of the Andreev transport 
for large values of U\. It implies that the effects of U\ 
is mainly to renormalize the inter-dot hopping matrix 
element t. In contrast, the Coulomb interaction U2 in 
the dot connected to the SC lead causes the qualitative 
changes, particularly near the crossover region between 
the Andreev behavior and Kondo behavior, as seen in 
Fig.EKa). 

We next study the precise feature of the crossover, 
choosing the Coulomb interaction at the QD1 to be 
U\ = and £1 = 0, for simplicity. In this case, the 
correlation effect of U2 on the low-energy properties can 
be described by a local Fermi liquid theory for a single 
impurity, and the fixed-point Hamiltonian for free quasi- 
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particles can be expressed in the form 



54 



jj(0) 



A rf2 (dLdl +H.c.) +t 



'£ (4 



CT d 2a + H.c. 



+ 6 (n d>2 - 1) + ^t.jv + i? 



(22) 



Here, t and Ad2 are the renormalized inter-dot coupling 
and the onsite SC potential in the QD2, respectively. The 
definitions of these renormalized parameters are provided 
in the appendix 151 Note that in a special case, for £2=0, 
the renormalized level position becomes £2 = 0, and then 
the conductance is determined by the ratio of t and Ts 
(= A d2 ), 



4e 2 



r N r s 



1 



r N r £ 



(23) 



The value of these renormalized parameters can be de- 
duced from the fixed point of NRGJ^ Particularly for 
£1 = £ 2 = 0, the ratio of the renormalized parameters 
appearing in Eq. (f2U|) links to the phase shift of Bogoli- 
ubov particles such that j?/(T N r s ) = —cotip. 

Figure [6] shows the NRG results of t, Tg, and the ra- 
tio ^/(Tn^s) as functions of i/IV for several values of 
{/2/rs. The parameter (t/t) 2 shown in Fig. H^a) corre- 
sponds to the wave function renormalization factor Z, 
which gives a measure of the correlation effect due to 
U2 (sec also the appendix [B]). For weak interactions 
U2/TS 1-5 the renormalization factor is almost con- 
stant (t/t) 2 ~ 1.0, and shows only a weak t/T^ depen- 
dence. It means that the energy scale of the resonance 
width is unrcnormalizcd. We have confirmed that these 
features are unchanged qualitatively in the parameter re- 
gion of U 2 /Ts < 2.0. The behavior changes significantly 
for U2/TS > 2.0. In this large U2 region, (t/t) 2 decreases 
with t/Tff, and goes to zero for t/T^ <C 1. This indicates 
that the system approaches a strongly correlated regime. 
In the opposite limit t/T^ 3> 1 the renormalization fac- 
tor gets close to (t/t) 2 — > 1, and for large t the ground 
state can be characterized by an inter-dot singlet state. 

The renormalized SC pair potential Ts (= Ad 2 ) is 
plotted in Fig. [(Jb). We see that Ts decreases as the 
Coulomb interaction [7 2 increases. Furthermore in the 
limit of t — > 0, Ts vanishes for strong interactions 
U2/TS > 2.0, while it reaches a finite value for weak 
interactions U2/TS < 2.0. Reflecting the behavior of t 
and Ts, the ratio 1? /(TnTs), which determines the con- 
ductance via Eq. (|2U)) . also shows a clear difference be- 
tween the strong and weak interaction regions. We see in 
Fig. He) that for U 2 /T s < 2.0 the ratio i 2 /^? s) de- 
creases monotonically as t/T^ decreases from 1.5 to 0. In 
contrast, ^/(TnTs) takes a minimum for U2/TS > 2.0 
at an intermediate value of t/T^, and then it increases 
as t decreases further. 




u 


2 /r s 







-0- 


1 




1.5 




2 




2.01 




2.5 


-A- 


3 




5 



t/T 



N 



FIG. 6: (Color online) Renormalized parameters (a) (t/t) 2 , 
(b) Ts/Ts and (c) the ratio of the parameters P/(TnTs) as 
a function of t/T at for several values of {/2/Ts. The other 
parameters are Ui = 0, £1 = £2 = and Tn = Ts- 
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FIG. 7: (Color online) Pair correlation at the QD2, («d,2) = 
{d\^d\^ + d 2 _id 2 ^), as a function of t for Ui = 0, £1 = £2 = 
and Tn = Ts- 



The difference between the weak and strong interaction 
regions can be seen also in the expectation values of the 
pair correlation at the QD2 defined by 

(K d>2 ) = (4t4| + d 2id 2t ) = I-Q2, (24) 

where Q 2 = (n 7j _2) is the number of the Bogoliubov 
particles in the QD2. Note that in the present case the 
Friedel sum rule described in Eq. (TT31 does not give the 
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local charge of a single site, so that we have calculated 
Q2 directly from the above definition. Figure [7] shows 
the results. In the weak interaction region U2 < 2Ts, 
the pair correlation (nd,2) in the QD2 increases as t gets 
smaller, and takes a value {nd,2) = 1 for t — > 0, which 
corresponds to the occupation number of the Bogoliubov 
particle Q2 = 0. In contrast, in the strong interaction 
region U2 > 2Ts, ( K d,2) has a maximum at finite t, and 
then it decreases to reach (red, 2) = for t —> 0, which 
corresponds to the occupation number Q2 = 1. 

The significant changes taking place at L^/Ts = 2.0 
can be explained as a result of the crossover in the ground 
state. In order to figure this crossover out, we consider an 
atomic limit t 0, where the QD2 is connected only to 
the SC lead. In this limit the eigenstates of the H de- 
fined in Eq. ([3]) can be obtained simply by diagonalizing 
a single site Hamiltonian of the QD2 with the onsite SC 
pair potential Ad2 = and the Coulomb interaction 
U2- Then, the ground state becomes a SC spin-singlet 
state for U% < 2Ts, and a magnetic doublet state for 
U2 > 2Ts- For small inter-dot coupling t, the essen- 
tial features seen in the atomic limit still remain espe- 
cially for the local SC singlet state for U2 < 2Ts- In 
the doublet region of U2 > 2Ts, however, the conduc- 
tion electrons, which come from the normal lead through 
the QD1, screen the local moment emerging in the QD2 
to form a Kondo singlet state. The level crossing which 
takes place at U2 = 2Ts in the atomic limit becomes the 
crossover for finite inter-dot coupling. 

These two singlet states correspond also to two differ- 
ent fixed points of NRG. The Hamiltonian H eS given in 
Eq. describes the interacting Bogoliubov particles, 
as mentioned. In the present case for £1 = £2 = 0, the 
energy level takes the value of E2 = T$- Therefore, the 
number of the Bogoliubov particles Q2 in the QD2 de- 
pends crucially on Ts and U2- The occupation number 
approaches the value of Q2 = for small t in the weak 
interaction region U2/2 < E2, and thus the QD2 becomes 
almost empty. This situation corresponds exactly to the 
frozen impurity fixed point, which is one of the basic fixed 
points for the asymmetric Anderson models Note that 
the frozen impurity fixed point for the Bogoliubov par- 
ticles corresponds to the local SC singlet state for the 
original electrons, as (nd.2) = 1 for Q2 = 0. On the other 
hand, in the strong interaction region of U2/2 > E2, 
the occupation number of the Bogoliubov particles ap- 
proaches the value of Q2 = 1 for small t — > 0. This 
situation corresponds to usual Kondo regime, which is 
referred to as the strong coupling fixed points 

From the above perspective, we can consider further 
the behavior of the conductance seen in Figs. [3ja), [4] 
and Efa). In the weak interaction region U2 < 2Ts, 
the Andreev reflection contributes to the tunneling cur- 
rent. This is because the renormalized SC pair poten- 
tial Ts (= A,j2) remains finite in this region, although 
the Coulomb interaction U2 suppresses the SC proxim- 
ity. The conductance takes the unitary- limit value 4e 2 //i 
on the resonance, which for £1 = £2 = takes place in the 
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FIG. 8: (Color online) Conductance as a function of the 
energy level £2 in the QD2 for si + Ui/2 = 0, Ui = U2 = 4T S 
and Fm = Ts- The arrows indicate the value of £2/^2 at 
which the singlet-doublet transition occurs in the atomic limit 
t -¥ 0. 

case that the ratio i /(rjvTg) becomes equal to 1. The 
resonance condition can be satisfied in the weak interac- 
tion region since T 5 and t are finite. The situation is quite 
different in the strong interaction region U2 > 2Ts, where 
the ground state is a Kondo singlet. For small inter-dot 
coupling t, the renormalized SC pair potential Ts is sup- 
pressed by the Coulomb interaction 1/2, and tends to zero 
for large L^- The renormalization factor Z (=t 2 /t 2 ) also 
approaches zero in the strong-coupling limit. Further- 
more, the ratio ^/^FaTs) has a local minimum in the 
strong coupling case as shown in Fig. [DJc). In a narrow 
parameter region near the crossover point, the conduc- 
tance has the two successive peaks (see the conductance 
for U2 = 2.0ir,5 in Fig. [S(a)). This is caused by the fact 
that the value of the ratio ^/(TnTs) at the local min- 
imum becomes less than 1. In the case of Fig. E|a), we 
have confirmed that the two-peak structure is seen in the 
region 2r s < U 2 < C/ 2 * with J7 2 * ^ 2.05r s . For U 2 > t/ 2 *, 
the conductance does not reach the unitary-limit value, 
as the minimum value of the ratio becomes greater than 
1. 



C. gate- voltage dependence 

Wc have seen in the above that the behavior of the con- 
ductance depends strongly on which interaction, U\ or 
U2, is changed. In this subsection, we discuss the effects 
of gate voltage on the Andreev transport. Particularly, 
we examine the role of the energy levels e\ in the QD1 
and £2 in the QD2 separately, choosing the interactions 
the same U\ = U2- 



1. £2 dependence 

Figure [5] shows the £2 dependence of the conduc- 
tance for several values of t/Ts- These results are ob- 
tained from Eq. Q16p , by choosing the energy level in the 
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FIG. 9: (Color online) Conductance as a function of the 
energy level e\ in the QD1 for e 2 + t/ 2 /2 = 0, Ui = U 2 = 4T S 
and Tjv = Vs- The symbols (*) on an upper horizontal axis 
correspond to the phase boundaries shown in Fig. 1101 




singlet 




-0.5 
Ej/Uj 



FIG. 10: (Color online) Phase diagram of the ground state 
of an isolated DQD system described by Hg S + Hdqd, for 
Ui = Ui = 4Fs and £2 + t/2/2 = 0. Symbols indicate the 
points of the single-doublet transition for t/Ts = 1.0, 1.5 and 
2.0. 



QD1 and the other parameters such that E\ = —Ui/2, 
Ui = U2 = 4rs and Tjy = Tg- We see that the conduc- 
tance has two peaks, which are broadened as the inter-dot 
coupling t increases. These two peaks reflect the level 
crossing which takes place at E2 = t/2/2 in the weak 
coupling limit t -> 0. Here, E 2 = y/{e 2 + U 2 /2) 2 + A} 2 
with Ad2 = Ts is the energy level of the Bogoliubov 
particles defined in Eq. (fTO]). The two arrows in Fig. [8] 
correspond to the critical points, and the ground state 
is a magnetic doublet in the region between the arrows 
—0.93 < £2/^2 ^ —0.07, while the ground state becomes 
a singlet on the outside. 

For finite inter-dot coupling t, the conduction elec- 
trons from the normal lead contribute to the screening 
of the local moment in the magnetic doublet state, and 
the ground state for —0.93 < £2/^2 < —0.07 becomes a 
Kondo singlet state.— Therefore, the conductance peaks 
seen in Fig. [5] are caused essentially by the crossover be- 
tween the Kondo singlet state and the local SC singlet 
state, particularly for small t. As t increases, the ground 
state varies gradually to an inter-dot singlet between the 
QD1 and QD2. 



2. £1 dependence 

We next consider the £\ dependence of the conduc- 
tance. To this end, we have to carry out the NRG cal- 
culations in the situation E\ + U±/2 ^ in the QD1 and 
A(j2 7^ in the QD2. In this case, the system no longer 
has the uniaxial pseudo-spin symmetry mentioned in Sec. 
Ill Bt and the number of Bogoliubov particles Q t ot is n °t 
conserved. For this reason, we have calculated the con- 
ductance using the Kubo formula given in Eq. (|17j) . 

Figure |9] shows the conductance at zero temperature 
as a function of £1 for several values of t/Ts, choosing 
the Coulomb interactions to be U\ = U2 = 4Ts and 
e 2 = —t/2/2. For weak inter-dot coupling of the value 
of t/Ts = 0.5, the ground state can be characterized 
by the Kondo singlet at —1.0 < £i/t/i < 0.0, and the 
conductance is suppressed by the Coulomb interaction. 



In contrast, for the strong coupling cases t/Ts > 1.0, the 
ground state near £\/U\ ~ —0.5 corresponds the inter- 
dot singlet as seen in Fig. [3] In this case, the conductance 
shows two peaks, and the peaks move away from the 
symmetric point at £\/U\ = —0.5 as t increases. In order 
to clarify these features, we have diagonalized a two-site 
Hamiltonian £t| ff + Hdqd consisting of the QD1 and 
QD2. The ground-state phase diagram of this cluster is 
plotted in Fig.[TO]in a £\/U\ vs t/Ts plane. The ground 
state is a spin-singlet in the central region between the 
two lines of the phase boundary, and is a spin-doublet 
state on the outside. Comparing this figure with Fig. |H1 
we see that the positions of the conductance peaks agree 
well with the phase boundary between the singlet and 
doublet ground states of the isolated DQD. 

We now consider the screening of the local moment 
emerging at the QD2 in the doublet regions seen in 
Fig. [TO] In these regions the energy level of the QD1 
is away from the Fermi level, and thus the QD1 is almost 
empty or doubly occupied. Therefore, the spin degree of 
freedom is frozen in the QD1, and it does not contribute 
to the screening. Nevertheless, the charge fluctuation still 
survives in the QD1, so that the conduction electrons 
from the normal lead can tunnel into the QD2 via the 
QD1 to screen the local moment. This kind of the virtual 
process induces an antiferromagnetic exchange coupling 
between the moment in the normal lead and QD2. This 
screening process is analogous also to a super-exchange 
mechanism, in which the conduction electrons come over 
a potential barrier. Therefore, the singlet bond which 
contributes to the Kondo screening becomes long in this 
case. This state has a deformed Kondo cloud, and is 
different from the singlet state that we discussed in Sec. 
IIII Al and IIII Bl because near half- filling the spin degree of 
freedom at the QD1 contributes to the screening. From 
these observations, it can be concluded that the conduc- 
tance peaks in Fig. [^correspond to the crossover between 
the inter-dot singlet and the long-range Kondo singlet. 
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IV. SUMMARY 

We have studied the Andreev transport through a 
DQD coupled in series to normal and SC leads, using the 
NRG method. We have demonstrated that in the case 
that the energy level of the QD adjacent to the normal 
lead is on the electron- hole symmetric point, namely for 
£i = 0, the system can be mapped onto a two-impurity 
Anderson model for the interacting Bogoliubov particles. 
In this case, the conductance at zero temperature is de- 
termined by the phase shift of the Bogoliubov particles. 

In order to clarify the ground-state properties, we 
have obtained the local SC correlation at the dots and 
the renormalizcd parameters for the low-energy quasi- 
particlcs. We have shown that the low-energy properties 
for a small t region are described by two distinct fixed 
points for the Bogoliubov particles. The local SC singlet 
state for small interaction U < 2T$ is described by the 
frozen-impurity fixed point, and the Kondo singlet state 
for U > 2Ts is characterized by the strong-coupling fixed 
point. Moreover, for large t, the two dots prefer to form 
the inter-dot singlet state. 

The conductance in the local SC singlet region for U < 
2Ts takes the unitary-limit value Ae 2 /h of the Andreev 
transport. In the Kondo region, we have found that there 
is a narrow region near the crossover point U — 2T$, in 
which the two-peak structure in the i-depcndcncc of the 
conductance appears. Outside this narrow region, the 
conductance peak in the Kondo regime gets smaller with 
increasing U. We have also found that for large t, the 
two dots form the inter-dot singlet, and it causes a low 
conductance in this regime. 

We have also calculated the gate- voltage dependence of 
the conductance, by varying the energy level in each dot 
separately. The conductance is enhanced at the transient 
region, where the crossover between two different types 
of the singlet states takes place. In particular, the energy 
level £1 of the QD1 causes the deformation of the Kondo 
screening cloud and makes the singlet bond long. Corre- 
spondingly, the conductance has the peaks which signify 
the crossover to this long-range Kondo screening. 

The results wc have obtained in the present work are 
for the large gap limit A — > 00. Therefore, the virtual 
processes using the continuum states outside of the SC 
gap A have not been taken into account. There are, gen- 
erally, some quantitative corrections for finite A. Nev- 
ertheless, since the perturbation expansion with respect 
to 1/A is applicable to the ground state properties, the 
results obtained in the large gap limit can be regarded as 
the zero-th order contributions. In fact, in a single QD 
coupled to a superconductor which corresponds to the 
limit t — > of the present system, the corrections due to 
finite A vary the ground-state phase diagram quantita- 
tively. Namely, the phase boundary of the singlet-doublet 
transition deviates from the one for the large gap limit 
A — > 00, and the spin-singlet region becomes wide for 
finite Ai 39 ' 55 i 56 A similar correction to the crossover en- 
ergy scale will also arise in our system. 



Recently, a DQD system coupled to superconductors 
has been fabricated using the carbon nanotubei^. The 
realization of such a system will provide further inter- 
esting examples of correlation effects in the context of 
the Andreev transport in nanoscale systems. We hope 
that the situation we have discussed in this paper will be 
examined experimentally in the near future. 
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Appendix A: Phase shift of the Bogoliubov particles 

We consider the Green's function for the Bogoliubov 
particles in the DQD, which is described by the Hamil- 
tonian defined in Eq. (TTTT) for £1 = 0, 



G_DQ_o(t — t') 



T 



1-2 At) 



Cj -j 



r-i(* 
r-i(* 



t ) G 7 _ 1j7 _ 



n 



(Al) 



Note that the matrix form of GoQD{t— t') stems not from 
the Nambu spinor but from a single component Green's 
function for the two interacting dots. The Fourier trans- 
formation of the retarded Green's function for Eq. (|A1[) , 
G r D Q D (e), can be represented as follows, 



{G r DQD (s)} *= 



-iT N -E r n {e) 



t-EJ 2 (e) 
e-E 2 -Z r 22 (e) 



(A2) 



where (e) is the self-energy due to the Coulomb inter- 
action. The advanced Green's function G a D Q D (e) is also 
given by taking the Hermitian conjugate of G r D Q D {e) 1 

namely, G a D Q D (e) = {G r D Q D (e)Y . Because of the time 
reversal symmetry of the Hamiltonian (fTTj) . £' 12 (£) = 
^21 i £ ) IS satisfied. At zero temperature, the imaginary 
part of the self-energy at the Fermi level e — van- 
ishes, i.e. ImS^(0) = 0. Therefore, the inverse matrix of 
G r D Q D (e) at e = can be expressed in the form, 



{G r DQD (0)}- = 



t* 



t* 



(A3) 



10 



where 



E 2 = E% 



^=HeE; i (0), 

t* = i-Re£[ 2 (0) 



ReS^ 2 (0), 



(A4) 



In the second equality, we have used the relation Eq. 
(US]). 

Appendix B: Renormalized parameters for Ui = 
and 6=0 



Using the Friedel sum rule for the Hamiltonian (|lip.— 
the local charge at the DQD for the Bogoliubov particles 
Q is given by 



2cp 

IT 



(A5) 



where (p is the phase shift of the Bogoliubov particles 
defined asi£ 



In the case the Coulomb interaction in the QD1 is zero 
U\ = with £1 = 0, it is possible to deduce the renormal- 
ized parameters for the Bogoliubov particles using NRG. 
This is because the self energies E^e) and Y7 l2 (e) be- 
come zero in this case, and the retarded Green's function 
G r D Q D (e) in eq. (|A2[) takes a simplified form 



det{G a DQD (0)Y 
det{G r DQD (0)}" 



7r — tan 



TnE* 2 



(A6) 



The local charge Q for the Bogoliubov particles is pro- 
portional to the pair correlation K and also to the local 
charge M for the electrons, as mentioned in Sec. Ill Bl 

Similarly, the conductance G in the case of £1 = can 
be expressed in terms of the phase shift <p. Using the 
Landauer formula, the linear conductance at zero tem- 
perature can be expressed in the form 2 ^ 



{G^(£)r = 



e + iTpf t 

t e-E 2 -V r 22 {e) 



(Bl) 



The asymptotic form of the Green's function at the QD2 
for small e ~ at zero temperature is given by 



Lt 7_2,7-2 



e-E 2 



t ■ 



e + iT N 



(B2) 



g = ^ 4 f 2 

h 



N 



F f- 



T,/- 



,|(0) 



(A7) where 



the Fourier transform of an 



Here ' Ff-tfj-^te) is 
anomalous Green's function for the electrons at the 
defined 



QD1, denned as 

0<{/-it(*)./-u(*')}} 

formation of Eq, 



■ F /_it,/_ 1 j.(* 



n 



f) = -iQ(t - 
Applying the Bogoliubov trans- 
I f i j_(0) can be rewritten as 



E 2 =Z [E 2 + E r 22 (0)\, t=y/Zt, (B3) 



e=0 



F 



/_it,/_i4. 



(0) = 2uvReG 



7-1.7- 



-x(0) 



(A8) 



Here, G r 1 _ l ~_ 1 (e) is the Green's function of the Bogoli- 
ubov particles at the QD1 defined in Eq. (|A1[) as the (1,1) 
element of G r D Q D (e), and the value at e = is given in 
Eq. (|A"3]l. Using this value with Eqs. © and (fAH) . the 
formula given in (|A7[) can be expressed in terms of the 
phase shift, 



G = 



4e 2 / A d n \ 2 4 \ 



±d2 

E 2 



( {t*Y-EUt \ 



1 + 



IT It 1 si " V 



The free quasi-particles, which are characterized by these 
renormalized parameters, particularly E 2 and t, can be 
described by an effective Hamiltonian which is given in 
Eq. (|22|) in terms of the electrons with the parameters 
defined by 



A d2 = 2uvE 2 = Z 



E5a(0) 



E-, 



A 



d2 



6 = (M 2 -f 2 )^ 2 = £ 



^ 2 (0) 



£ 2 



(A9) Here, we have used u and v defined in Eq. ©. 



(B5) 
(B6) 
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